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Abstract 

Enhanced proximity effect by mid gap Andreev resonant state (MARS) in a diffusive normal 
metal / insulator / triplet superconductor (DN/TS) junction is studied based on the Keldysh- 
Nambu quasiclassical Green's function formalism. By choosing a p- wave superconductor as a typical 
example of TS, conductance of the junction and the spatial variation of quasiparticle local density 
of states (LDOS) in DN are calculated as the function of the magnitudes of the resistance Ra, 
Thouless energy in DN and the transparency of the insulating barrier. The resulting conductance 
spectrum has a zero bias conductance peak (ZBCP) and LDOS has a zero energy peak (ZEP) 
except for a = 7r/2 (0 < a < 7r/2), where a denotes the angle between the lobe direction of the 
p-wave pair potential to the normal to the interface. The widths of the ZBCP and ZEP are reduced 
with the increase of Rd while their heights are drastically enhanced. These peaks are revealed to 
be suppressed by applying a magnetic field. When the magnitude of Rd/Ro is sufficiently large, 
the total zero voltage resistance of the junction is almost independent of the Rd for a 7^ 7r/2. 
The extreme case is o; = 0, where total zero voltage resistance is always Rq/2. We also studied 
the charge transport in + ipy-w&ve junctions, where only the quasiparticle with perpendicular 
injection feel MARS. Even in this case, the resulting LDOS in DN has a ZEP. Thus, the existence 
of the ZEP in LDOS of DN region is a remarkable feature for DN/TS junctions which have never 
been expected for DN/ singlet superconductor junctions where MARS and proximity effect compete 
with each other. Based on these results, a crucial test to identify triplet pairing superconductors 
based on tunneling experiments is proposed. 
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I. INTRODUCTION 



Nowadays, enormous numbers of unconventional superconductors have been discovered 
where pair potential changes sign on the Fermi surface [l|. 




1^. It is known that re- 
tunneling spectroscopy of 
1^. The most dramatic ef- 



flecting on the sign change of the pair potential 
unconventional superconductor is essentially phase sensitive 
feet is the appearance of zero bias conductance peak (ZBCP) in normal metal / insulator 
/ unconventional superconductor junctions due to the formation of mid gap Andreev res- 
onant states (MARS) lOl 11 1. The ZBCP are observed in the actual many experiments 
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The MARS is a unique resonant state expected for t 
perconductor assisted by the Andreev reflection 0, 10, 



le interface of unconventional su- 



28|. The origin of this MARS is 



due to the anomalous interference effect of quasiparticles at the interface where injected 
and reflected quasiparticles from the unconventional superconductor side feel different sign 



of the pair potentials 



port properties 
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It is well known t 
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lat MARS influences several charge trans- 
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Since MARS is expected both for triplet and singlet superconductor junctions j48|, |49|, it 
is a challen ging is sue to present a new idea to discriminate triplet superconductin g st ates 

00,0,0,13,0 

from singlet ones through tunneling spectroscopy via MARS 



Although there are several studies about unconventional superconductor junctions, 

^ 3, 0, 0, 3, 0, 0, 0, Q 0, H Q H 0, H S 
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80| almost all of them are restricted to ballistic regime. If we take into account impurity 



scattering, we can expect various interesting features even for conventional superconductor 



junctions jSJJ 
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93| |. As regards diffusive normal 



metal / insulator / conventional singlet s-wave superconductorfDN/ CSS) iunctions, there 



has been a reniarkable progress in theories of proxim i ty e f fect 



mi\M 



P 94. 



effect in unconventional superconductors |ll4 1 has been yet to be clarified. Recently, we have 
developed a theory of charge transport in diff usive norm al metal / insulator / unconventional 

exte nding Na zarov's theory of matrix 



01 proxirn i ty e i tect WJdL 

i3, Ii3, Iml, Him 
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However, proximity 
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singlet s upercond uctor (DN/USS) junctions 
current |lOll Ill2| within quasiclassical treatrnent 
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applied to ci-wave superconductor junctions |118| . Unfortunately, however, it is revealed 
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that the proximity effect and the MARS compete with each other in DN/USS junctions. 
Although the interface resistance is reduced by the MARS irrespective of the magnitude of 
the transparency at the interface, the resuhing total resistance of the junction R is always 
larger than Rq/2 + R^, where Rd and Rq is the resistance in DN and Sharvin resistance at 
the interface, respectively. This is because that the angular average of many channels at 
the DN/USS interface destruct the phase coherence of the MARS and the proximity effect 



[see Fig. 1. of ref. 



5J]). This destructive angular average is due to the sign change of 



the pair potentials felt by quasiparticles with injection angle and those with —0, where 
the angle is measured from the direction normal to the junction interface. However, 
in diffusive normal metal / insulator / triplet superconductor (DN/TS) junctions, we can 
escape from the above destructive average. In the last paper, we have presented a theory 
which is available for DN/TS junctions 5^. It is revealed that charge transport in diffusive 
normal metal / triplet superconductor (DN/TS) junctions is significantly unusual. We can 
expect enhanced proximity effect by the MARS. The total zero voltage resistance R in the 
DN/TS junctions is significantly reduced by the enhanced proximity effect in the presence 
of the MARS. It is remarkable that when the resistance in DN, Rd, is sufficiently larger 
than the Sharvin resistance Rq, R is given hj R = Rq/C^, ©which can become much smaller 
than the preexisting@lower limit value of R, i.e., Rq/2 + Rd. In the above, C_ is a constant 
completely independent of both Rd and Ri,, where Ri, denotes the interface resistance in 
the normal state. When all quasiparticles injected at the interface feel the MARS, R is 
reduced to be i? = -Ro/2 irrespective of the magnitude of Rd and Rh. At the same time, 
local density of states (LDOS) in the DN region has zero energy peak (ZEP) due to the 
penetration of the MARS into the DN region from the triplet superconductor (TS) side 
of the DN/TS interface. These outstanding features have never been expected either in 
DN/CSS junctions or DN/USS junctions. However, the reference js^ does not contain 
the necessary technical details of the matrix current derivation and the obtained results 
are limited. In order to understand this remarkable enhanced proximity effect by MARS 
peculiar to the DN/TS junctions much more in detail, one has to evaluate the conductance 
spectrum and quasiparticle density of states in wide range of several parameters. 

In this paper, we present a detailed derivation of the matrix current in (DN/TS) junctions. 
We express compact formula of matrix current and Keldysh-Nambu(KN) Green's functions 
relevant to the actual boundary condition in comparison with those of DN/USS junctions. 
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Here, we restrict our attention to triplet superconductors with Sz — 0, where Sz denotes the 
z component of the total spin of a Cooper pair. In order to show up the anomalous charge 
transport in DN/TS junctions due to the coexistence of the MARS and the proximity effect, 
we present detailed numerical calculations of the conductance spectra of DN/TS junctions 
for p-wave superconductors. We investigate the dependence of the spectra of the bias voltage 
conductance on various parameters: the height of the barrier at the interface, resistance Rd 
in DN, the Thouless energy Eth in DN and the angle a between the normal to the interface 
and the lobe direction of p-wave superconductor. Eth can be expressed by Eth — TiD/ LF' 
with diffusion constant D in DN and the length of DN. We normalize the voltage-dependent 
conductance as{eV) by its value in the normal state, a^^ so that ax{eV) = as{eV)/(Xi^. At 
the same time, wc focus on the spatial dependence of LDOS in DN for DN/TS junctions. We 
also studied about several related physical quantities which identify the anomalous charge 
transport. Our main results are as follows: 

1. The ZBCP is always seen in the shape of arieV) for a ^ 7r/2 with < a < 7r/2. The 
magnitude of (Tr(0) is drastically enhanced with the increase of Rd/Rb- The half width of 
the ZBCP, which is proportional to exp(— Cc-R(i/-Ro) for C'c, is almost independent of the 
transparency of the DN/TS interface and E^h- On the other hand, aT{eV) for DN/USS 
junctions with dxy-^s^re superconductor shows a very different behavior. The magnitude of 
arieV) is reduced with the increase of Rd/Rb due to the absence of the proximity effect. 

2. The LDOS in DN has a ZEP except for the case with a = 7r/2 where MARS is absent. 
The existence of the ZEP in DN means that the penetration of the MARS into DN. The 
height of the ZEP is significantly enhanced with the increase of Rd/ Rq- The half width of 
the ZEP is proportional to exp{—CpRd/ Rq) where Cp is a constant almost independent of 
Z. 

3. We can express the LDOS in the DN region using the proximity parameter 9 as Real[cos 9]. 
If we denote the 9 at DN/TS interface as ^o, = 2iRdCOsa/Ro is satisfied at £ = where 
quasiparticle energy e is measured from the Fermi energy. Since is a pure imaginary 
number for £ = 0, we can expect a ZEP in LDOS. This unique feature has never been 
expected for DN/USS of DN/CSS cases, where £ = is always a real number. 

4. The total zero voltage resistance R in the DN/TS junctions is significantly reduced by 
the enhanced proximity effect in the presence of the MARS. It is remarkable that when Rd is 
sufficiently larger than the Sharvin resistance Rq, Ris reduced to be = -Ro/(2 cos a) , which 
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can become much smaller than the preexisting lower limit value of i?, i.e., Rq/2 + Rd. For 
low transparent junctions, R is also reduced to be i? = Ro/{2 cos a). When all quasiparticles 
injected at the interface feel the MARS, R is reduced to be i? = i?o/2 irrespective of the 
magnitude of Rd and i?^. This drastic situation is realized for a = 0. 

6. The sharp ZBCP or ZEP in LDOS due to the enhanced proximity effect is sensitive to the 
applied magnetic field H. The height of ZBCP is significantly reduced by H. The threshold 
value of the magnetic field is Hxh — 8h/{eSDN) where Sdn denotes the magnitude of the 
area of DN region. 

7. As a candidate of TS (triplet superconductor), we also choose px + ipy-wave supercon- 
ductor. This superconducting state is so called chiral superconducting state where broken 
time reversal symmetry state (BTRSS) is realized. Although only quasiparticles with per- 
pendicular injection feel MARS, both LDOS and (7T{eV) has ZEP and ZBCP due to the 
enhanced proximity effect by MARS. 

It should be remarked that these remarkable features have never been expected either 
in DN/CSS or DN/USS junctions. The structure of the paper is as follows. We formulate 
the model in use in section 2. We also present there the detailed derivation of the matrix 
current and end up with the expression for the normalized conductance. We focus on p-wave 
superconductor junctions in section 3 and evaluate arieV), p{e), zero voltage resistance and 
the measure of the proximity effect 6q for various cases. The last part of this section, we 
discuss Px + ipy-wave junctions. We summarize the results in section 4. 



II. FORMULATION 



In this section we introduce the model and the formalism. We consider a junction consist- 
ing of normal and superconducting reservoirs connected by a quasi-one-dimensional diffusive 
conductor (DN) with a length L much larger than the mean free path as in our previous 
paper. The interface between the DN conductor and the TS (triplet superconductor) has a 
resistance Rb while the DN/N interface has zero resistance. A schematic illustration of the 
model we use are shown in Fig. ^ The positions of the DN/N interface and the DN/TS 
inter face are denoted as a; = —L and x = 0, respectively. According to the circuit theory 
|ll2j |. the constriction area {—Li < x < Li) between DN and TS is considered as com- 
posed of the diffusive isotropization zone {—Li < x < —L2), the left side ballistic zone 
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(— ^2 < X < 0), the right side balhstic zone (0 < x < Li) and the scattering zone (x = 0). 
The scattering zone is modeled as an insulating delta-function barrier with the transparency 
T = T(0) = 4 cos^ 0/(4 cos^ + Z"^), where Z is a dimensionless constant, is the injec- 
tion angle measured from the interface normal to the junction. Here, we express insulating 
barrier as a delta function model Hi,6{x), where Z is given hj Z = 2mHi,/ (h'^kp) with 
Fermi momentum kp and effective mass m. In order to clarify charge transport in DN/TS 
junctions, we must obtain Keldysh-Nambu (KN) Green's function, which has indices of 
transport channels and the direction of motion along x axis taking into account the proper 
boundary conditions. We restrict our attention to triplet superconductors with Sz = that 
preserves time reversal symmetry. Sz denotes the z component of the total spin of a Cooper 
pair. It is by no means easy to formulate a charge transport of DN/TS junctions since the 
quasiparticle Green's function has no angular d ependence by the impurity scattering in the 
DN. Ho wever , as shown in our previous paper [iItI llisl - if we concentrate on the matrix 
, 112! via the TS to or from the DN, we can make a boundary condition of the 



currents 



KN Green's function. The sizes of the ballistic zone in the DN and t 
the current flow direction are much shorter than the coherence length 
we assume the fiat interface, the mome ntum para llel to the interface is conserved. Thus it is 



le s c attering zone in 
11811. Since 



112 



possible to construct a matrix current 
in TS. 
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112| based on the asymptotic Green's function 
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FIG. 1: Schematic illustration of the model of constriction area. 
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In this section, we will show how to derive the matrix current in DN/TS junctions. Then 
we will derive the retarded and Keldysh components of the matrix current. Finally, we will 
show how to calculate the conductance of the junctions. 



A. Usadel equation in DN 

In order to clarify the charge transport in DN/TS junctions, we first concentrate on the 
KN Green's function in DN. We define KN Green's function in DN as 

/ 

Gn{x) = 



Rn{x) Kn{x) 
iiv(^) 



(1) 



Here, we are using quasiclassical approximation. Then, we can express Rn{x) as 

Rn{x) — cos t/j sin 9ti + sin ip sin 9t2 + cos 9ts, (2) 

with Pauli matrix in the electron hole space, fi, T2, and 73. Since Rn{x) obeys Usadel 
equation, following equations are satisfied, 

hDl^e - i^f cos e sin 9] + 2ie sin ^ = 0, (3) 
ax^ ax 

In the present case, since no supercurrent is flowing in the junctions, |^ = is satisfled. 
We can choose ip = tJjQ where iIjq is a constant independent of x. The boundary condition of 
Gn{x) at DN/TS interface is given by 

L dGNjx) _ < 7(0) > 

/(0)=2[Gi,5] 

B = (-TifGi, HZ^] + HZ^H+ - T^GiHZ^H+Gi)-\Ti{l - HZ^) + T^GiHZ^H+) (5) 

with Gi = Gi^{x = 0_), H± = {G2+ ± G2-)/2, and Ti = T/{2 -T + where G2± is 

the asymptotic Green's function in TS as defined in our previous papers. In the above, the 
average over the various angles of injected particles at the interface is defined as 

< I {(f)) ci0 cos 0/(0)/ ci0r(0)cos0 (6) 

J—k/2 J-n/2 



with = I and T{(f)) = T. The resistance of the interface Rb is given by 

= (7) 

I-ir/2 d<PT{4>) COS0 

with Sharvin resistance at the interface, Rq. In the above, Gi and G2± can be given by 




(8) 



R2± K2± 

\ isi 

where the Keldysh component Ki,2± is given by Ki(2±) = -Ri(2±)/i(2)(0) - /i(2)(0)Ai(2±) 
with the retarded component -Ri,2± and the advanced component Ai^2± using distribution 
function /i(2)(0). In the above, i?2± is expressed by 

R2± = {.g±rz + f±f2) 



with g± = e/ye?^ — A^, f± = A±/ sj A% — e"^, and A2± = — f3_R2±'r3 where e denotes the 
quasiparticle energy measured from the Fermi energy. A+ (A_) is the effective pair potential 
felt by quasiparticles with an injection angle (tt — 0) as shown in Fig. [TJ /2(0) = /o5(0) = 
tanh[£:/(2kBT)] in thermal equilibrium with temperature T. Here, we put the electrical 
potential zero in the TS. We also denote H^, H^, B, I = I{(f)) as follows. 



f R ( 



Y 




H-= , (9) 





B 



B. Calculation of the retarded part of the matrix current 

First, we pay attention to the boundary condition of the retarded part of KN Green's 
function at DN/TS interface in order to determine the value of ipo. The left side of the 
boundary condition of eq. Q can be expressed by 

—Rn{x)^Rn{x)\x=o = — [-sin?/'ofi + cosijjof2]{—)\x=o (10) 

rid Ox rid OX 

In the above, the most remarkable point is that the fs component is vanishing. In order to 
calculate the right side of eq. (jH), we must study the parity of Ir as a function of 0. In 



general, can be expressed by using several spectral vectors, 



Ir = AiTiidR ■ dR)-'{--{l + T^){S2+ - S2-)^[si X (S2+ + S2-)] ■ T 



+2TiSi ■ (S2+ X S2-)[Sl X (S2+ X S2-)] " T 
+2riSi ■ (S2+ - S2-)[Sl X (S2+ - S2-)] " t 
-i{l+T^){l - S2+ • S2_)[Si X (S2+ X S2-)] " T 
+ 2iTi(l - S2+ ■ S2_)[Si ■ (S2+ - S2_)Si - (S2+ - S2-~ 



;i + T2)(s2+ X S2-) - 2TiSi X (S2+ - S2-) " 2T^Si ■ (S2+ X S2-)Si 



(11) 

(12) 



with i?i = Si • T and R2± = S2± ■ t. 

The spectral vectors si and S2± are given by 



^ sin 6^0 cos 



Si 



sin 6'o sin ^/'o 





f ° 1 




/±(0) 







(13) 



where denotes the 6* at x = 0_. We postulate that component of < Ir > should be 
zero. Here, we focus on the parity of Ir since there is an angular average over (p in the actual 
calculation. As a comparison, we also look at the case where TS (triplet superconductor) 
is substituted with USS (unconventional singlet superconductor). For DN/TS junctions, 
9±{~(p) = Qt^'P) and f±{—4>) = — /=p(0) are satisfied, while for DN/USS junctions, g±{—(p) = 
(7zp(0) and f±{—4>) = /=f(0) are satisfied. {cIr ■ (Ir)'^ can be given by 



[dR ■ cIr) 



^ \ rfW (0) + rfW (0) sin iPo DN/TS 



1 



rfW(0) 



DN/USS 



(14) 



for DN/TS and DN/USS junctions, where 4'\4>) and dW(0) are even functions with 6, 
while d^^^ ((f)) is an odd function with 0. We apply similar discussions for other terms in Eq. 



II Bl). For the convenience, we define 







(15) 
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After simple manipulations, we can show the following relations 



^(1 + T^){S2+ - S2-)'[Sl X (S2+ + S2-)] " S, 

Fi^?(0)cos^o DN/TS 
Fil\(l)) cos ipo DN/USS 



(16) 



2riSi • (S2+ X S2-)[Si X (S2+ X S2-)] • S, 

(0) COS Vo sin V^o DN/TS 
F^e^ (0) cos V'o sin V'o DN/USS 



(17) 



2T1S1 • (S2+ - S2-)[Sl X (S2+ - S2-)] ' 

(Fi*)(0)sin7Ao + Fi?)cosV^o DN/TS 
f£^(0)cosV'o DN/USS 

- i(l + Ti^)(1 - S2+ • S2-)[Sl X (S2+ X S2-)] " 

Fi*)((/))sinV'o DN/TS 
F£\(/.)sin^/^o DN/USS 



(18) 



(19) 



2iTx{l - S2+ ■ S2-)[Si ■ {S2+ - S2-)Si - {S2+ - S2-)] ■ S, 



(20) 



Fi?(0)sin^o + i^i?(</') DN/TS 
F5^^\0)sin^o DN/USS 

In the above, Fi?(0), ^i? ^2?(</') and F^'Xcf^), are even func- 

tions of 0, while FiJ{(j)), F^^J{(j)), ^4*^(0), F^^J{(j)), and -^50^(0) are odd functions of 0, 
respectively. By postulating that the fs component of Ir should vanish after the angular av- 
erage over 4>, we can show that siu'^o = for DN/TS junctions and cos^^o = for DN/USS 
junctions. 

The resulting retarded part of the boundary condition is given by 

< F > 



L 09 



Rh 



'^T{fs 00800 - gs sin^o) 
2-r + r(cos^o5s + sin^o/s) 



(21) 
(22) 
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with gs = + + + fs = i{Ug--g+f-)/{l + g+g- + Uf-) for DN/TS 

junctions and gs = (^+ + g~)/{l + ^+^- + /+/-), fs = (/+ + /-)/(! + 9+9- + /+/-) for 
DN/USS junctions, where g± = g±{4>) and f± = f±{4>). This is one of the central results of 
this paper. 

C. Calculation of the Keldysh part of the matrix current 

Next, we focus on the Keldysh component. We define as 

h = ^Tr[f3/x] Ik = 2{R,Bk + KiBa - BrK, - BkA,) (23) 

with K, = rJ,{0) - A(0)ii /i(0) = /ojv(O) + /3^(0)f3. Br is given by 

b?n + hfh + hfh DN/TS 
h^^h + l)i^T2 + fei'Va DN/USS 



B 



R 



,(t) Ti{Tismeo + fs) -(t) Tifs (<) Ti{Ti cos Oo + gs) 

— 1 , Oo — , 0\' — T , 

1 A ' ^ A ' A 

,(s) TJs As) 7i(Tisin6'o + /s) As) Ti{Ti cos Oq + gs) 

Ol = 7 , On — -. , Oo = . 

1 A ^ A A 

A = (1 + T^) + 2Ti(^5 cos ^0 + fs sin ^o) 

In the above gs, fs and fs are defined by 



9s = < 



{g+ + g.)/{l + g+g. + f+f.) DN/TS 
(^+ + ^_)/(l + 5+^_ + /+/_) DN/USS 



fs={ 



i{f+g^ - f-g+)/{l + g+g- + /+/_) DN/TS 
(/+ + /_)/(! + g^g^ + /+/_) DN/USS 



fs^{ 



(24) 



(25) 



(26) 



(27) 



(/+ + /_)/(! + g+g. + /+/_) DN/TS 
i(/+^_ - g+f-)/{l + g+g- + /+/_) DN/USS 

We can calculate Ba similar to the case of Br. After simple manipulation, we can show 

that Ba = —f^Bl^fs and the resulting is given by 

h = Trace{f3(^i^i^ + ^I^i^) (28) 
-[hiRlBl + BrR^ + + R\BR)]foj,{Q) - [(4 + R{){Br + si)]/3Ar(0)}/2, 

12 



After angular average over (f) and the operation of the trace, it is easily shown that the 
second term which is proportional to /oAr(O) does not contribute to < /(, >. It is necessary 
to obtain Bk which is given by 

Bk = Dr^Nk - DrDkBa (29) 

with Ba ^ Da^Na and 



D = -Ti[Gi, HZ^] + HI^H+ - T^GiHZ^H+Gi, D 

where Nk and Na is the Keldysh and advanced part of N given by 

N = Ti- TiHZ^ + TlGiHZ^H+ N = 




Dr Dk 
Da 



We can express Nk and Dk as linear combinations of distribution functions fos{^), /oAr(O), 
and /3jv(0) as follows both for DN/TS and DN/USS junctions. 



N^ 



K 



C?fos{0) + C?foN{0) + (7£V3iv(0), DN/TS 
cSfosiO) + cSfoNiO) + cSfMO), DN/USS 



(30) 



^ . ci'^fosiO) + C'?Voiv(0) + (7?V3Ar(0), DN/TS ^^^^ 
cSfosiO) + C'£Voiv(0) + (7^?/3iv(0), DN/USS 

, eg + DN/TS ^ ^ 

^ii = -^r/ (32) 

C}^^ DN/USS 

Wltn — L/j^g -|- L^ig , L/2 — L/2e "T L/2o , '_/4 — L/4g -|- , L/5 — Ly^^ -f- L/50 , dUQ 

(7^ ^ (AW ^ (j(t)^ hj 2x2 matrix C]^^ with r = t, s, j = 1, ..,7 and k = e, o. When the 
suffix /c is e (o), C*]^"* is an even (odd) function of 0. , , C^] , C^^] , Cgg , and C^g 
are hnear combinations of i and f2, while , C20 ^ , C^] , , and C*7g are linear 
combinations of fi and fa. 1 is a unit matrix in the electron hole space. On the other hand, 
Ci^J , C20 and are linear combinations of f2 and f3, while C'^o 1 C'io , C'^^o 1 and C'fo are 
linear combinations of 1 and fi. Taking account of these facts, after angular average over 0, 
< lb > can be given by 

< lb > = < 'IVaCe{(^i + ^I)[^;^Ef3 - {Br + Si)]/3Ar(0)} > /2, 

13 



Bke = D],'[C3-C,D2'Na], 

- T,{R,r, - fgiOi^^ipii - T,R,R;^'Rp{R,rs - fail)]. (33) 

with Cs = C^} ((73 = Ci'^) and Cq = C^^ {C^ = C^'J) for DN/TS (DN/USS) junctions. 
Since following equations are satisfied, 



(34) 



Dji^iTi — TiR^l + T^RiR^Rp) — I3r, 



(35) 



-B_R,(1 + Rjn^ + TiRiRpRj;) — TiR^Rp 



(36) 



< Jb > is given by 



<h>=^< Trace{-(i?i + R\)Br{Ri + - (^i + Ri){Br + Bj^)} > fsN{0) 

For the later convenience, we define < Iho > with < I^o >=< h > //3Ar(0). Then the final 
resulting expression of < ho > is given by the following equation, 

T Co 



< ho >=< 



2 I {2-T)+T{coseogs + smeofs) P 



> 



(37) 



Co = T(l+ I cos^o r + I sin 00 \')[\ gs \' + \ fs I' +1+ I fs 
+4(2 - T)[Real(^5)Real(cos0o) + Real ( ) Real ( sin ^o)] 
+4T[Imag(cos 9o sin eQ)lmag{fsg*s)] ■ 



In the above, the definition of gs, fs and fs are given in eqs. 
6o = 0, we can reproduce well known results in ballistic limit 




and (P7j) . By choosing 
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D. Calculation of the conductance 



The electric current is expressed using as 



'"-im-J,, ''^^\MG.(.)^n (38) 

where (Giv(x)^^^)^ denotes the Keldysh component of (Giv(x)^^^). In the actual 
calculation, we introduce a parameter 6 = 6{x) which is a measure of the proximity effect 
in DN as described in the previous subsections, where we denoted 9{0) = Oq. Using ^(x), 
Rn{x) can be denoted as 

Rn{x) = f^cos6{x) + f2^\TL0{x). (39) 

An{x) and Kn{.x) satisfy the following equations, Ai\]{x) = —f3R\^{x)f3, and Kn{x) = 
f^Ar(x)/i(x) — fi{x)Ai\i{x) with the distribution function which is given by = 

/oAr(x)+f3/3jv(a;). In the above, fsN^x) is the relevant distribution function which determines 
the conductance of the junction we are now concentrating on. From the retarded or advanced 
component of the Usadel equation, the spatial dependence of 9{x) is determined by the 
following equation 

hD—e{x) + lie sin[^(x)] = 0, (40) 
ox^ 

while for the Keldysh component we obtain 

with 6imag{x) = lmag[0(x)]. At X = —L, since DN is attached to the normal electrode, 
9{—L)=0 and fs^^—L) = fto are satisfied with 

fto = i{tanh[(£ + eV)/{2kBT)] - tanh[(£ - eV)/{2kBT)]}, 

where V is the applied bias voltage. As shown in our previous paper, from the Keldysh part 
of Eq. (jU, we obtain 

|-(^)cosh^Imag(^o) U^o_= (42) 
After a simple manipulation, we can obtain /37v(0_) 

RbftO 



/3iv(0- 



p, _|_ Rd<ho> fO dx 
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Since the electric current lei can be expressed via in the following form 




we obtain the following final result for the current 




(43) 



Then the total resistance R at zero temperature is given by 



< IbO> L J-L COSh^ 9imag{x) 




(44) 



In the following section, we will discuss the normalized conductance arieV) — 
as{eV)/aN{eV) where o's(N){eV) is the voltage-dependent conductance in the superconduct- 
ing (normal) state given by as{eV) = 1/R and aN{eV) = = ^/{Rd + Rb), respectively. 

It should be remarked that in the present theory, R^/ R^ can be varied independently 
of T, i.e. of Z, since we can change the magnitude of the constriction area independently. 
In the other words, Rd/Rb is no more proportional to Tav{L/l), where Tav is the averaged 
transmissivity and I is the mean free path in the diffusive region, respectively. Based on this 
fact, we can choose Rd/Rb and Z as independent parameters. 

III. RESULTS 

In this section, we focus on the line shapes of the conductance of DN/TS junctions where 
p-wave symmetry is chosen as a pairing symmetry of triplet superconductor (TS). The pair 
potentials A± are given by A± = ±Ao cos[(0 ^ a)] where a denotes the angle between the 
normal to the interface and the lobe direction of the p-wave pair potential and Aq is the 
maximum amplitude of the pair potential. In the above, denotes the injection angle of 
the quasiparticle measured from the x-axis. In the following, we choose < a < 7r/2. It 
is known that quasiparticles with injection angle with — 7r/2 -|-Q;<0<7r/2 — o; feel the 
MARS at the interface and induce ZBCP. 

A. Line shapes of the conductance 

Let us first choose a = where all quasiparticles feel MARS. For alpha = 0, we also call 
Pj;-wave case in the following. In this case, at eV — 0, the total resistance of the junction 
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Ris Ro/2 and completely independent of Rd as shown in our previous paper j5J| . Thus the 
resulting ariO) is o"t(0) = 2{Rd + Rb)/ Ro- In Figs. |21and|ni voltage dependent conductance 
axieV) is plotted for low (Fig. |2)) and high transparent (Fig. ^ interface. For Z = 3 (Fig. 
12)), we always expect ZBCP independent of the detailed value of Rd/Rb- With the increase 
of the magnitude of Rd/Rb, the width of the ZBCP becomes narrow (see curves b and c in 
the left panel), while it's height is drastically enhanced. The width of the ZBCP is rather 
reduced by choosing small magnitude of i?Th(see curves b and c in the right panels of Fig. 

EI). 



15i ' ' ' I ■ ■ ■ 1 ■ ■ ■ I ■ ■ 
10- -- 




-0.1 0.1 -0.1 0.1 

eV/Ao eV/Ao 



FIG. 2: Normalized conductance aT{eV) for Z=3, and a = 0. Eth = for the left panel and 
Erh = O.O2A0 for the right panel, respectively, a, Rd/Rb = 0; b, Rd/Rb = 0.1; and c, Rd/Rb = 1- 

Similar plots for Z = is shown in Fig. |21 In this case, Rb = Rq is satisfied due to 
the absence of the insulating barrier at the interface. For Rd/Rb = 0, aT{eV) has a broad 
ZBCP due to the absence of the normal reflection at the interface P]. With the increase of 
the magnitude of Rd/Rb, the width of the ZBCP is reduced (see curves b and c both in the 
left and right panels), while its height is enhanced. As seen from the right panel of Fig. |2l 
the width of the ZBCP is reduced drastically with the decrease of the magnitude of Ej^h- 

In Fig. EJ the corresponding plots for DN/USS (diffusive normal metal / insulator / 
unconventional singlet superconductor )junctions are shown with rf-wave superconductor, 
where the misorientation angle between the normal to the interface and the crystal axis 
is chosen to be /3 = 7r/4. For (3 = it /A, we also call d^^-wave case in the following. A± 
is given by A± = ±Aosin(2^). Although quasiparticles feel MARS independent of their 
injection angles, as shown in our previous papers 11^1 . the proximity effect and MARS 
strongly compete with each other and the proximity effect is completely absent. aT{eV) is 
independent of the value of Eth- Both for Z = and Z = 3, the magnitude of axieV) is 
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reduced monotonically with the increase of Rd/Rb- This feature is completely different from 
that in DN/TS junctions as shown in Figs. |21and|ni 

Next, we look at a dependence of axieV) for DN/TS junctions. In this case, an injected 
quasiparticle with the injection angle (p with —njl + a<(f)<'K/2 — a feel the MARS. It is 
shown in our previous paper, that only the quasiparticle channel with this injection angle 
can contribute to the proximity effect. Only at a = 7r/2, neither proximity effect nor MARS 
exist. As shown in Fig. ZBCP is absent for tt/2. In other cases, aT{eV) has a ZBCP. 

In order to clarify the proximity effect in DN, it is necessary to focus on the local density 
of states (LDOS) of the quasiparticles in DN region. In Fig. IHl LDOS of DN/TS junctions 
is plotted for a = with Eth = and Eth = 0.02Ao. 

In Fig. [7| the corresponding plot with various a for Ej-h = Aq is shown. p{e) is the 
normalized LDOS by its value in the normal state, where e denotes the quasiparticle energy 
measured from the Fermi surface. The curves a, b and c denote the LDOS at x = —L/A, 
—L/2 and —L, respectively. Since DN is connected to the normal electrode at x = —L, 
p{e) = 1 is satisfied independent of e as shown in curves c in Figs. IHlandlZI As shown in Fig. 
El the p{e) has a zero energy peak (ZEP) in DN (curves a and b in left and right panels). 
Even if a deviates from 0, ZEP in LDOS does not vanish (curves a and b in the left and 
middle panels of Fig. [7j). This is because that quasiparticles with injection angle with 
—7i/2 + a < (f) < 7r/2 — a feel MARS and can contribute to the proximity effect. Exceptional 
case is a = 7r/2, where proximity effect is absent. Then the resulting p{e) = 1 independent 
of e. 

In order to compare with DN/USS cases, we choose (i-wave pair potential with A± = 

4 



-0.4 0.4 -0.4 0.4 

eV/Ao eV/Ao 

FIG. 3: Normalized conductance axieV) for Z=0, and a = 0. E^h = for the left panel and 
Erh = O.O2A0 for the right panel, respectively, a, Rd/Rb = 0; b, Rd/Rb = 0.1; and c, Rd/Rb = 1- 
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FIG. 4: Normalized conductance axieV) for DN/USS junctions with d^y-wave superconductor. 
Z = 3 for the left panel and Z = for the right panel, a, Rd/Rb = 0; b, Rd/Rb = 0.1; and c, 
Rd/Rb = 1. 



Aocos[2(^^ =1= f3)], where f3 denotes th e ang le between the normal to the interface and the 



crystal axis of d-wave superconductor 



118[. The corresponding p{e) is shown in Fig. |S| p{e) 



is nearly constant as a function of e and does not have a sharp ZEP (Figs. ElandlTj). In 
this case, although quasiparticles with injection angle tt/A — P <\ |< 7r/4 + /? feel MARS, 
they can not contribute to the proximity effect. For = 0, MARS is absent and proximity 
effect becomes conventional one (see also Fig. 1 of |ll8|). p{€) at x = — L/4 has a gap like 
structure (curve a in the left panel of Fig. |H1). Although p{e) at x = — L/4 has a broad 
peak like structure for (3 = vr/S, p(0) < 1 is satisfied contrary to the DN/TS junction's 
case. For (3 = n/A, although aT{eV) has a ZBCP, due to the absence of the proximity effect, 
p{e) = 1 for any case. Thus we can conclude that line shapes of p{e) in DN region of DN/TS 
junctions are significantly different from those of DN/USS junctions. 
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-0.4 0.4 -0.4 0.4 
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FIG. 5: Normalized conductance aT{eV) for Z=3, and Rd/Rb = 1- E^h = for the left panel 
and Exh = 0.02Ao for the right panel, respectively, a, a = 0; b, a = '7r/4; and c, a = tt/2. 



19 




FIG. 6: Normalized local density of states p{e) in DN for Z=3, Rd/Rb = 1 and a = 0. Exh = 
(left panel) and Eth = 0.02Ao (right panel), a, x = — L/4; h, x = —L/2; and c, x = —L. 
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FIG. 7: Normalized local density of states p{e) in DN for Z=3, Rd/Rb = 1 and ii^j'/j = Aq. a = vr/4 
(left panel), a = Svr/S (middle panel), and a = 7r/2 (right panel), a, x = — L/4; b, x = —L/2] and 
c, X = — L. 

In Figs, iniandlini the proximity parameter Oq = 9{x = 0_) is plotted as a function of e. 
In Fig. HU a dependence of the Oq is plotted. The magnitude of Real(Imag)(6'o) is increased 
with the increase of Rd/Rb- Contrary to the case of DN/USS junctions, Real(6'o) vanishes 
at e = while Imag(6'o) remains nonzero at e = 0. For fully transparent case (Fig. EJ, 
Real(6'o) and Imag(6'o) decrease with the increase of e/Aq for Eth = Aq. For Exh = 0.02Ao, 
Real(6'o) decreases from and has a minimum at about e ~ 0.02Ao. Imag(6'o) decreases 
rather rapidly as compared to the case with Eth = Aq. For low transparent case with Z = 3 
(Fig. ITn|l . Real(^o) has a sudden change around e ^ 0. This sudden change is remarkable for 
the large magnitude of Rd/Rb (see curves b and c). With the increase of the magnitude of 
a, the resulting magnitude of Oq is reduced as shown in Fig. [TT] For a = 7r/2, the resulting 
9o is zero. 

As seen from the sudden change of near zero energy for low transparent junction, we 
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FIG. 8: Normalized local density of states p{e) in DN for DN/USS junction with Z = 3, Rd/Rb = 1 
and Exh = ^o- We have chosen d-wave superconductor with A-t = Aocos[2(0 =p (3)]. /3 = (left 
panel), (3 = vr/S (middle panel), and (3 = 7r/4 (right panel), a, x = — h, x = —L/2; and c, 
X = —L. 
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FIG. 9: Real(upper panels) and imaginary parts(lower panels) of ^ire plotted as a function 
of e for a = and Z = with E't^/j = Aq (left panels) and E^h = 0.02Ao (right panels), a, 
i^d/i^b = 0.1; b, Rd/Rb = 0.5; and c, Rd/Rb = 1- 

can imagine that there is a very small energy scale. Actually, as seen from Fig. El the 
width of axieV) is significantly reduced for low transparent junction with the increase of 
the magnitude of Rd- The sudden reduction of the magnitude of axieV) around eV ~ is 
related to the sudden change of 9q around e = 0. 

In order to understand this small energy scale, we focus on the half width of the zero 
bias conductance peak of arieV). We define the half width Ec as ariEc) = ^o"t(0) and 
p{Ep) = |p(0), respectively. In Fig. [T21 Ec is plotted as a function of Rd/ Rq for various Z. 
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FIG. 10: Real(uppcr panels) and imaginary parts(lower panels) of 9q arc plotted as a function of 
e for a = and Z = 3. Exh = (left panels) and Eth = 0.02Ao (right panels), a, Rd/Rb = 0.1; 
b, Rd/Rb = 0.5; and c, Rd/Rb = 1- 
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FIG. 11: Real(uppcr panels) and imaginary parts(lower panels) of are plotted as a function of 
£ for Rd/Rb = 1 and Z = 3. Eth = •^o (left panels) and Eth = 0.02Ao (right panels), a, a = 0; 
b, a = 0.257r; and c, a = O.STStt. 

If the magnitude of Rd/Ro is large, Eq can be expressed by Ec/A.o ~ Eqi exp{—CcRd/Ro): 
where Eqi and Cc are independent of Rd- It is interesting that the magnitude of Cc is 
almost independent oi Z. As seen from left and right panels, the magnitude of Cc is almost 



22 



constant by the change of Exh- It is a very unique feature of the present proximity effect that 
Ec has an exponential dependence of Rd and the magnitude of Ec is drastically reduced 
with the increase of Rd- The corresponding p{e) also has a sharp ZEP as shown in Fig. IHI 
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FIG. 12: Ec is plotted as a function of Rd/Rb for = Aq (left panel) and Exh = 0.02Ao (right 
panel), a, Z = 1; b, Z = 2; and c, Z = 3. 

In order to understand this small energy scale, we define Ep as p{Ep) = |p(0). In Fig. [T31 
Ep is plotted as a function of Rd/ Rq- Ep can be expressed by Ep/ Aq ~ Epi exp{—CpRd/ Rq) 
for the large magnitude of Rd/ Rq where Epi and Cp are independent of Rd- Comparing the 
curvature of curves from a to c, the magnitude of Cp is almost independent of Z- As seen 
from left and right panels, the magnitude of Cp is almost constant by the change of Eth- 
The magnitude of Ep is drastically suppressed with the increase of Rd as in the case of Ec- 
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FIG. 13: Ep is plotted as a function of Rd/Rb for E^h = (left panel) and Ej^h = 0.02Ao (right 
panel), a, Z = 1; b, Z = 2; and c, Z = 3. 



B. Properties at zero voltage 

In this subsection, we focus on the zero voltage properties of DN/TS junctions with p- 
wave superconductor, where we can obtain several analytical results without solving Usadel 
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equation numerically. The Oq is given by by solving eq. dH). By using eqs. ()37|) and 

the zero voltage resistance R is given by 

^ ^ ^ tanh^ 2 

°^ /i [1 + exp(2 I 9o^ \)]h + cosh2^oJ2 ^ ^ 

/7r/2— a 
COS 
-n/2+a 

-7r/2+a (0) COS /■'^/2 2T2 (0) COS 



-./2 [2-T(0)]2 - y./2-« [2-T(0)]2 - 
with ^^oi = 2RdCosa/Ro. First, we calculate the zero- voltage resistance {e 0) at different 
values of a as a function of Rt/Ro for the DN/TS junctions with two extreme cases (curves 
a and b of Fig. IT^ . For a = 7r/2, namely the p^-wave case, R increases linearly as a function 
of Rd, where no proximity effect appears (curve a of Fig. HH). For a = 0, Px-^ave case, 
where R = Ro/2 is satisfied independent of Rd (curve b of Fig. ITlj) . This anomalous R 
dependence is a most striking feature by the enhanced proximity effect by the MARS. The 
corresponding result for the DN/CSS junctions with s-wave (curve c) and DN/USS junctions 
with dxy-wave case (curve d) is also plotted as a reference. For s-wave case, it is well known 
jo^ that R has a reentrant behavior dR/dRd |ij^=o< as shown in curve c of Fig. (fT^ . In 
p-wave cases, this reentrant behavior of R does not appear. Actually, from the relation 

OR 

" Rd=0 



OR 



D 



cos^ a 



1 : 2. (46) 

cos a -\- J^/2^a [2-T(</.)]2 

we can show that dR/dRo |_r^=o> is always satisfied. For dxy-wave case, due to the 
formation of the MARS for all as in the case of pa;-wave junction, R/Rb = Ro/{2Rh) at 
Rd = 0, which is identical to that for p^-wave junction (curve d of Fig. IT^ . However, for 
nonzero Rd, R/Rb = Ro/{2Rb) + Rd/ Rb is satisfie d, due to the absence of the proximity 
effect as discussed in our previous papers IiItI lll^l . 

The a dependence of R is shown in Fig. ^1 For Z = 0, as shown in curves b and c, R 
increases gradually with Rd/Rb- However, for the large magnitude of Rd/Rb, the ratio of 
the increment is saturated and R is almost constant as a function of Rd/Rb- For Z = 3, i? 
is nearly constant with the increase of Rd/Rb- As seen from Eq. ()45|) . for sufficiently large 
magnitude of Rd/ Rq, the second term in Eq. becomes negligibly small and R converges 
to -Ro/(2cosa) independent of the magnitudes of Z and Rb- On the other hand, for small 
transparent junctions, i.e., the large magnitudes of Z, we can neglect the term proportional 
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FIG. 14: Total zero voltage resistance of the junctions R is plotted as a function of Rd/Rb for 
Z = 1.5 with a, a = 0.257r; and b, a = 0.47r. The curve c and d presents the same dependence 
for the DN/CSS junctions with s-wave superconductor and DN/USS junctions with d^y-wave 
superconductor, respectively. 

to I2 in Eq. (jlSj) . In such a case, R converges to be i?o/(2cosa) independent of Rd- For 
Z = 3 (see the right panel of Fig. IT^ . R saturates to be constant much more rapidly with 
the increase of Rd/Rb as compared to the corresponding case for Z = (see the right panel 

of Fig. nni). 
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FIG. 15: Total zero voltage resistance of the junctions R is plotted as a function of Rd/Rb for 
Z = (left panel) and Z = 3 (right panel), a, a = 0; b, a = 7r/2; and c, a = OAn. 

Before closing this subsection, we pay attention to the zero voltage electrostatic potential 
distribution at zero temperature (poi^)- Due to the enhanced proximity effect by the MARS, 
the spatial dependence of 4>o{x) is very unusual. The electrostatic potential in the DN region 
is expressed by 

^i^) = 1^ deTTiice\KN{x)] = / (i£/3Ar(a;)Real[cos 6']. (47) 
4 JO Jo 
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The zero-bias voltage static potential distribution is defined by [10 



Since /3Ar(x) is expressed by 



D _|_ Rd<IbO> fO dx 

r / \ _ L Jx cosh'^ e.magix) , . 

hN{X) - Ra<I,Q> rO dx /ffl' i'^yj 

with 9imag{x) = lma.g[9 (x)] , the resulting (f)o{x) is given by 

_j_ Rd<ho> fO rfa; 

0o(x) = Real[cos^^(x)]--^^^|^-;^^^ (50) 

at zero temperature. In Fig. UHl x dependence of (Pq{x) is plotted both DN/USS junction 

with drcy-wave superconductor and DN/TS junction with p^.-wave superconductor. Although 

quasiparticles always feel MARS both for dxy-v/ave and p^^-wave cases, for dxy-v/ave case the 

proximity effect is completely absent, while for ^^.-wave case, proximity effect is enhanced 

by MARS. For d^y-^ave case, since 9{x) = 0, (poi^) is given by 

LRq — 2xRd , . 

= LW + 2fiJ- 

Actually, as shown in Fig. (f)o{x) has a linear x dependence and decreases with x. This 
linear dependence corresponds to the ohm's rule. The magnitude of (f)o{x) is much reduced 
with the increase of Rd/Rb (see curves c and d). For pa;-wave case, since 9{x) = 2iRd/ Rq is 
satisfied, the resulting 4>o{x) is 

0o(x)=exp( — ). (52) 

L Kq 

Due to the enhanced proximity effect by MARS, 0o(^) is reduced drastically with the increase 
of Rd/Rb (see curves c and d). 

C. Magnetic field dependence of the conductance 

From the experimental view point, it is interesting to clarify how ZBCP in DN/TS junc- 
tions are influenced by the applied magnetic field H. It is known for conventional supercon- 
ductor junctions, that coherence of electrons by the prox imity effect is broken by the applied 
magnetic field or the magnetic impurity in DN 0, Il02| as shown in the following equation 



o2 ^ 

hD—9 + 2iesm9 sin 2^ = (53) 

OX^ Th 
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FIG. 16: Zero voltage electrostatic potential distribution at zero temperature (j)o{x) is plotted as a 
function of x for Z = 3. DN/USS junction with dxy-wave superconductor (left panel) and DN/TS 
junction with pa;-wave superconductor (left panel), a: Rd/Rb = 0, b: Rd/Rb = 0.1, c: Rd/Rb = 1; 
and d: Rd/Rb = 2. 



with h/rn = Qw^D^H^, where w is the width of the DN region and H is the applied magnetic 
field. With the increase of the magnetic field, the height of the ZBCP is suppressed as shown 
in Fig. El The sharp peak structure around the zero voltage is significantly suppressed 
when the magnitude of Ji/th is almost IOE't/i- This is because the magnitude of Im^o is 
significantly reduced by the increase of the magnitude of h/rn. 
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FIG. 17: Normalized conductance aT{eV) for DN/TS junctions with Z = 1.5 and Rd/Rb = 1 
where p-wave superconductor with a = is chosen. We choose Eth = for the left panel and 
Erh = O.lAo for the right panel, respectively, a, h/rH = 0; b, H/th = Eth', and c, h/rn = lOEnph- 

In order to understand the suppression of the proximity effect, we also look at p{e) in 
the DN region. The height of the ZEPs in p{e) both in the left [Ej-h = Aq) and right 
panels {Eth = 0.1 Aq) is reduced seriously with the increase of the magnitude of h/rn- For 
Ti/th = IOEth, pi^) — 1 for both cases (see curves c in the left and right panels in Fig. 
[TH|) . In this case, the proximity effect is almost absent and the ZBCP in axieV) in curves 
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c in Fig. El originates purely from the formation of the MARS at the DN/TS interface. 
The corresponding situation is reahzed in DN/USS junction with dxy-^scve superconductor 
(see Fig. E}. Here, we estimate the threshold magnetic field which suppresses the sharp 
ZBCP or ZEP originating from the enhanced proximity effect. Here, Ti/th and Eth are 
given by e^w^DH"^ / (Qh) and hD/L^, respectively. We can define threshold magnetic field as 
lOExh = e'^w'^DH^f^J (Qh). We can easily show that Hxh — By choosing L = Ix 10"^[m] 
w = 1 X 10~^[m], the resulting Hxh is 5 x 10"^ [Tesla]. 



Z=1.5 

ETh=Ao ETh=0.1Ao 




-0.05 0.05-0.1 0.1 



FIG. 18: Normalized local density of states of quasiparticle p{eV) for DN/TS junctions with 
Z = 1.5 and Ra/Rb = 1 where p-wave superconductor with a = is chosen. We choose x = — L/4. 
Eth = for the left panel and Exh = 0.1 Aq for the right panel, respectively, a, H/th = 0; b, 
H/th = Eth; and c, H/th = WExh- 



D. Conductance for p+ip-wave superconductor 

In this subsection, we discuss the case where Px + ipy-wave superconductor is chosen 
as TS(triplet superconductor). Stimulated by intensive studies about Sr2Ru04 after its 
discovery S], + ipy-wave superconductor, which is a p-wave superconductor with broken 
time symmetry, has a wide interest at the present. In this case, one of the important feature 
is that the MARS is formed only by the quasiparticle with perpendicular injection. The 
position of the resonance energy strongly depends on the injection angle of the quasiparticles. 
Also in this case, we can calculate the corresponding conductance axieV) following similar 
calculations as in section 2. Since the derivation of the matrix current and the relevant 
< Ibo > is rather long, we will o nly show the final results. The details of the calculations 
will be shown in elsewhere jl2c| . First, we look at the case with = 0, i.e., ballistic 
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junctions. As shown in curves a in the left and right panels of Fig. (TT(eV^) has a rather 
broad peak . This is because that the position of the resonance energy deviates with the 
increase of the magnitude of the quasiparticle energy e. For Exh = ^O; the magnitude of 
(Jt{^) is suppressed with the increase of Rd/Rb, contrary to the case with DN/TS junctions 
without BTRSS as shown in Figs. |21 El and El For E^h = 0.02Ao, with the increase of 
Rd/Rb, o"t(0) has a sharp and narrow peak around zero voltage. This feature becomes 
prominent with the increase of Rd/Rb- 
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FIG. 19: Normalized conductance axieV) for DN/TS junctions with px + ipy superconductor for 
.Z^ = 3 is plotted. We choose Eth = for the left panel and Exh = 0.02Ao for the right panel, 
respectively, a, Rd/Rb = 0; b, Rd/Rb = 0.1; and c, Rd/Rb = 1- 

The corresponding plot for p{e) is shown in Fig. (201 p(^) is the normalized LDOS by its 
value in the normal state, where e denotes the quasiparticle energy measured from the Fermi 
surface. The curves a, b and c denote the LDOS at x = —L/4, x = —L/2 and x = —L, 
respectively. Since DN is connected to the normal electrode at x = —L, p{e) = 1 is satisfied 
independent of e as shown in curves c in the left and right panels of Fig. (201 p{e) has 
a zero energy peak (ZEP) in DN (curves a and h). Although only the quasiparticles with 
perpendicular injection at the DN/TS interface feel the MARS, the ZEP in LDOS remains 
as in the case in Figs. (HlandlZl The line shapes of LDOS for p^^ + ipj^-wave case is significantly 
different from the corresponding ones for DN/USS junctions with rf-wave superconductor 
(see Fig. The existence of the ZEP in LDOS of DN is a remarkable feature peculiar to 
triplet superconductor. We think that this peculiar property which has never been expected 
in singlet junctions may serve as a guide to identify the triplet superconducting state. 
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IV. CONCLUSIONS 



In the present paper, detailed theoretical investigation of the voltage-dependent conduc- 
tance of diffusive normal metal / insulator / triplet superconductor (DN/TS) junctions is 
presented. We have provided the detailed derivation of the formula of the matrix current 
presented in our previous paper For the reader's convenience, we explicitly present 

the retarded and the Keldysh parts of the matrix current in the DN/TS (diffusive normal 
metal / triplet superconductor ) junctions. Applying our general formula to DN/ insulator 
/ p-wave superconductor junctions, we have obtained the following main results. In the 
present paper, by changing the barrier parameter at the interface Z, the resistance Rf, at the 
DN/TS interface, the resistance Rd in DN, the Thouless energy Ejnh in DN and the angle 
a between the normal to the interface and the lobe direction of p-wave superconductor, we 
have studied the charge transport in DN/TS junctions in detail. 

1. The zero bias conductance peak (ZBCP) is always seen in the line shape oiaT{eV) except 
for a = 7t/2. The magnitude of ctt{0) is drastically enhanced with the increase of Rd/Rb- 
On the other hand, aT{eV) for diffusive normal metal / insulator / unconventional singlet 
superconductor (DN/USS) junctions with d^y-wave superconductor shows a very different 
behavior. The magnitude of arieV) is reduced with the increase of Rd/Rb due to the absence 
of the proximity effect. For the large magnitude of Rd/Ro, the half width of the ZBCP is 
proportional to exp{—Cc Rd/Ro), where Rq is the Sharvin resistance at the DN/TS interface. 
Cc is almost constant and independent of Z and Eth- 
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FIG. 20: Normalized local density of states p{e) for DN/TS junctions with + ipy-wave super- 
conductor is plotted for Z = 3. We choose Exh = for the left panel and E^h = 0.02Ao for the 
right panel, respectively, a, x = —L/A; b, x = —L/2; and c, x = —L. 



30 



2. The LDOS has a ZEP except for the case with a = ti/2 where MARS is absent. The height 
of the ZEP is significantly enhanced with the increase of Rd/Ro- For the large magnitude of 
Rd/Ro, the half width of the ZEP is proportional to exp{—CpRd/Ro)- Cp is almost constant 
and independent of Z and Eth- 

3. The proximity parameter at the DN/TS interface 9q is 2i cos aRd/Ro aX e = where 
quasiparticle energy e is measured from the Fermi energy. This unique feature has never 
been expected for DN/USS or diffusive normal metal / insulator / conventional s-wave 
singlet superconductor (DN/CSS) case, where £ = is always a real number. 

4. The total zero voltage resistance R in the DN/TS junctions is significantly reduced by 
the enhanced proximity effect in the presence of the MARS. It is remarkable that when Rd 
is sufficiently larger than the Sharvin resistance Rq, R is reduced to be i? = i?o/(2cosa), 
which can become much smaller than the preexisting lower limit value of R, i.e., Ro/2 + Rd- 
For low transparent junctions, R is also reduced to be i? = i?o/(2cosa) for any Rd- When 
all quasiparticles injected at the interface feel the MARS, R is reduced to be i? = Ro/2 
irrespective of the magnitude of Rd and Rb. This dramatic situation is realized for a = 0. 

5. The sharp ZBCP or ZEP in LDOS due to the enhanced proximity effect is sensitive to the 
applied magnetic field H. The height of ZBCP or ZEP is significantly reduced by H. The 
threshold value of the magnetic field is Hxh — 9>li/{eSDN) where Sdn denotes the magnitude 
of the area of DN region. 

6. We also choose px + ipy-wave state as a model of chiral superconductor. Although only 
quasiparticles with perpendicular injection feel MARS, crr(eV^) has a ZBCP and LDOS has 
a ZEP, due to the proximity effect. 

These novel features have never been expected either in DN/CSS or DN/USS junctions. 
Here, we propose an experimental setup which discriminate triplet superconducting state 
from singlet one by Scanning tunneling spectroscopy (STS) as shown in Fig. |^ Only 

/ /- / / A / 

FIG. 21: Schematic illustration of the experimental setup to discriminate triplet superconducting 
state from singlet one. LDOS measured by STS only have a ZEP for DN/TS junctions. 

for triplet case, LDOS measured by STS has a ZEP. The spatial dependence of the height 
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of ZEP can be observed by the proposed STS. According to our calculation discussed in 
the last section, we can obtain the LDOS for various x, where the positions of the DN/TS 
interface is denoted as x = 0. As in our previous section, we assume that DN is attached 
to the normal electrode (N) at x = —L. The spatial dependence of p{6) for e = is plotted 
in Fig. 1221 both for p^-wave and Px + ipy-wave cases. Both for two cases, the magnitude 




FIG. 22: Spatial dependence of the local density of state in DN region at zero energy is plotted 
for pa;-wave (left panel) and px + ipy-wave (right panel) superconductors with Z = 1. For p^^-wave 
case (left panel), Rd/Rb = 1 (curve a), and Rd/Rb = 0.5 (curve b). For + ipy-wave case (left 
panel), Rd/Rb = 2 (curve o), and Rd/Rb = 1 (curve b). Since we are concentrating on e = case, 
p{0) is independent of Et^- 

of p(0) decreases monotonically with the increase of the magnitude of the absolute value of 
x/L, where the position of the tip of STS changes from the DN/TS interface to the DN/N 
interface. We hope above spatial dependence of LDOS will be observed near future. As 
far as we know, the experiment of Knight shift measuring the uniform susceptibility is an 
only promising way to detect the triplet superconducting state. Our proposal to identify the 
triplet superconducting state based on the proximity effect in the presence of the MARS is 
an innovational idea. 

In the present paper, we have considered the normal metal / triplet or singlet supercon- 
ductor junctions, where the mean free path in the normal metal region is much shorter than 
the length of it. Thus, it is possible to use Usadel equation. In the actual junctions, there is 
a possibility that the magnitude of the mean free path in normal metal becomes as the same 
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order as its length. To describe this situation, we have to solve Eilenberger equation where 
the strengt h of t he impurity scattering can be changed as an input parameter. Very recently, 



Lofwander 



12ll | calculated the local density of states in normal metal / rf-wave superconduc- 



tor junctions based on the Eilenberger equation. His obtained result is consistent with our 
theory when the impurity scattering effect becomes significant. It is an interesting problem 
to study the charge transport in triplet junctions based on the Eilenberger equation. 

It is also possible to calculate tunneling conductance by the numerical simulation in 
the lattice model, where impurity potential is introduced in DN, where conductance is 



calcu^ 
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ated for vario us samples which have different configuration of randomness each other 



123 



124 . 



125| . Although it is difficult to take the large size of the system, the merit of 



the simulation is to take into account the impurity scattering effect exactly. The prominent 
property, ie., R is co mple tely independent of Rd for px-wave junction, is verified by the 



numerical simulation 



125|. 



A direct evidence of the mesoscopic interference effect by the proximity effect in uncon- 
ventional superconductors has been reported in size dependence effect of high- Tc supercon- 
ductor junctions. When the junction size becomes smaller, the conductance spectra are 
shown to be modified from the ballistic features. Since the high sensitivity of the ZBCP 
to the appli ed field was observed, the experimental data are consistent with the present 



results 



126[. Another way to detect the phase coherence in DN is to observe the LDOS 



of normal conductor in the vicinity of unconventional superconductors by tunneling spec- 
troscopy. There have been presented two different types of tunneli ng re sults on the LDOS 



observation of Au coated on YBCO. Experimental results of Ref. |127l | show the presence 
of gap structure both for Au/YBCO(100) and Au/YBCO(110), while Asulin et al. observed 
ZBCP on Au/YBCO(110) when Au i s thi nner than 7nm, while the ZBCP suddenly disap- 



pears as the Au lay er be comes thicker 128l |. These experimental reports are consistent with 
the present theory llisl ] in the point that MARS can not penetrate into DN in the case 
of ci-wave superconductors. However, since both of these measurements were performed on 
planer bilayer system (not mesoscopic size junctions), and the Au overlayer cannot be re- 
garded as DN, the detailed comparison of these results with the present theory seems to be 
inadequate. On the other hand, it is a challenging issue to make a junction using Sr2Ru04 
where ZBCP is already reported in tunneling spectroscopy [19]. In order to observe the 
enhanced proximity effect, we think the measurements in extremely small residual magnetic 
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field environments are necessary. This is because the ZBCP is easily suppressed even by the 
terrestrial magnetism. 

There are several problems which have not been discussed in the present paper. In the 



present study, we have focused on N/S junctions. The extension o: 
diffusive S/N/S junctions has been performed by Bezuglyi et al. 



the circuit theory to long 



lll|. In S/N/S junctions, 



the mechanism of multip 



on I-V curves 
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e Andreev ref l ections pr oduces the subharmonic gap structures 
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132, M3 



134 Il35| and the situation becomes much more 



complex as compared to N/S junctions. Moreover, in S/N/S junctions with unconventional 
superconductors, MARS leads to the anomalous current-phase relation and temperature 
dependence of the Josephson current [37]. An interesting problem is an extension of the 
circuit theory to S/N/S junctions with unconventional superconductors. Recent numerical 
simulation indicates the existence of the anomalous current phase relation in triplet junctions 



13q . In the present paper, since we follow the quasiclassical Green's function formalism, 
the impurity scattering is taken into account within the self- consistent Born approximation. 
It is a challenging problem to study the weak localization effects. 
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